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ABSTRACT 


The  question  whether  the  measurement  of  radar  cross-sections  at 
different  frequencies  provides  a  useful  diagnostic  tool  for  ascertaining 
the  electron  density  distribution  of  spherically  symmetric  plasma 
clouds  is  investigated.  This  is  accomplished  by  comparing  the 
calculated  radar  cross-sections  of  characteristic  plasmas  with  increas¬ 
ing  and  decreasing  refractive  index.  Exact  analytical  expressions  for 
the  radar  cross-sections  of  several  typical  plasma  spheres  with  in¬ 
creasing  and  decreasing  refractive  index  are  calculated.  The  calcu¬ 
lations  are  based  on  an  exact  wave  treatment  of  the  scattering  problem. 
Part  I  of  the  present  report  contains  the  exact  analytical  results 
obtained,  while  Part  II  will  be  devoted  to  the  numerical  evaluation  of 
these  results,  as  well  as  to  asymptotic  expressions  for  the  limiting 
cases  of  high  and  low  frequencies. 
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I. 


INTRODUCTION 


In  many  areas  of  practical  importance,  it  is  of  great  interest  to  be  able 
to  ascertain  the  electron  density  distribution  of  inhomogeneous  plasmas 
which  are  not  accessible  to  direct  measurement.  For  example,  such 
plasmas  or  clouds  of  ionization  might  arise  from  artificial  perturbations 
in  the  ionosphere  or  from  nuclear  detonations  in  the  atmosphere.  An 
accurate  knowledge  of  the  electron  density  distribution  in  such  plasmas 
would  be  of  great  aid  in  understanding  the  phenomena,  e.g.  the  nature 
of  the  reactions  and  their  respective  reaction  rates,  taking  place  inside 
the  plasmas. 

The  ultimate  purpose  of  our  present  investigation  is  to  determine  to 
what  extent  information  concerning  the  electron  density  distribution  of 
plasmas  can  be  obtained  by  means  of  ground-based  radar  cross-section 
measurements  carried  out  at  different  frequencies.  To  begin  with,  we 
shall  be  concerned  only  with  inhomogeneous  plasma  spheres  with  a 
spherically  symmetric  electron  density  distribution.  Accordingly,  we 


necessarily  limit  our  considerations  to  non-turbulent  plasmas.  More¬ 
over,  vve  shall  assume  that  all  electromagnetic  quantities  (e.g.  dielectric 
constant,  conductivity,  refractive  index)  describing  the  plasma  are 
likewise  spherically  symmetric  scalars.  We  are  thus  restricting  our¬ 
selves  to  the  case  where  the  background  magnetic  field  of  the  Earth 
with  its  attendant  anisotropy  is  negligible.  From  a  practical  point  of 
view,  our  results  will  then  be  valid  for  plasmas  for  which  the  electron 
collision  frequency  greatly  exceeds  the  Larmor  frequency  within  the 
plasma. 

Thus,  the  general  problem  to  which  we  address  ourselves  is  the  deter¬ 
mination  of  the  radar  cross-sections  of  inhomogeneous  (albeit  spherically 
symmetric)  plasmas  embedded  in  a  medium  of  uniform  electromagnetic 
properties.  This  outside  medium  need  not  necessarily  possess  the 
electromagnetic  properties  of  the  vacuum;  our  treatment  is  equally 
applicable  to  a  plasma  cloud  situated  in  the  ionosphere,  as  long  as  the 
electromagnetic  properties  of  the  surrounding  ionosphere  can  be  con¬ 
sidered  as  approximately  uniform. 

Our  approach  to  the  problem  will  consist  of  an  exact  wave  treatment 
of  the  scattering  of  a  plane  wave  from  the  plasma  under  consideration, 
based  on  Maxwell's  equations  without  the  introduction  of  any  approxima¬ 
tions  . 
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While  the  overall  motivation  of  the  present  research  is  thus  to  ascertain 
whether  the  measurement  of  radar  cross-sections  can  be  considered  a 
useful  diagnostic  tool  for  determining  the  electron  density  profile  of 
spherically  symmetric  plasmas,  the  present  study  does  not  attempt  to 
answer  this  question  in  its  most  general  form.  We  shall  be  concerned 
here  only  with  a  preliminary  investigation  designed  to  yield  a  comparison 
cf  the  radar  cross-sections  of  spherical  plasmas  with  increasing  and 
decreasing  complex  refractive  index,  as  a  function  of  radial  distance  from 
the  origin.  The  results  of  this  investigation  will  show  whether  the  two 
characteristic  cases  of  increasing  and  decreasing  refractive  index  give 
rise  to  distinctly  different  radar  cross-sections  as  a  function  of  fre¬ 
quency,  such  that  we  may  conclude  that  the  measurement  of  radar  cross- 
sections  holds  promise  as  a  diagnostic  tool  and  therefore  warrants 
further  investigation.  Future  planned  work  includes  the  possibility  of 
increasing  the  information  obtained  about  the  electron  density  profile  of 
plasmas  by  means  of  measurements  carried  out  by  ground-based  receivers 
at  different  locations.  This  would  involve  a  calculation  of  the  full 
differential  scattering  cross-sections  of  the  plasmas.  In  the  present 
preliminary  work,  however,  we  shall  be  concerned  only  with  proper 
radar  cross-sections. 
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In  Figure  1,  we  have  plotted  schematically  the  refractive  index  n  as 
a  function  of  radial  distance  r  for  the  four  general  cases  we  shall 
consider.  A  few  explanatory  remarks  may  be  called  for.  To  begin 
with,  the  plots  in  Fig.  1  are  only  schematic,  since  n  is  in  general  a 
complex  quantity  (we  include  absorption  in  our  analysis).  In  all  cases, 
the  space  outside  the  sphere  (r  >  b)  is  considered  to  have  uniform 
electromagnetic  properties  (n  =  constant).  Cases  (a)  and  (b)  represent 
the  extreme  cases  where  the  electrons  are  concentrated  with  constant 
density  in  the  center  and  in  the  outside  of  the  spherical  plasma,  res¬ 
pectively.  Cases  (c)  and  (d)  represent  a  more  continuous  variation. 

In  both  cases  (c)  and  (d),  the  presence  of  the  discontinuity  in  n  at 

r=a  represents  the  most  general  case  considered.  The  analysis  carried 

out  below  includes  the  special  cases  of  no  discontinuity;  thus,  we  have 

also  included  the  cases  where  in  case  (c),  n(a)  =  n  ,  and  in  case  (d), 

o 

n  decreases  continuously  to  the  constant  value  of  the  outside  medium 
at  r=0  ,  as  shown  by  the  dotted  curve.  Finally,  while  we  have  shown 
the  refractive  index  as  greater  than  unity  in  the  schematic  diagrams  of 
Fig.  1,  our  analysis  applies  equally  well  to  refractive  indices  which 
are  less  than  unity  or  negative. 

The  subject  matter  of  our  investigation  conveniently  separates  into  two 
portions.  Parti,  constituting  the  present  report,  is  concerned  with 


the  derivation  of  analytical  expressions  for  the  radar  cross-sections 
of  the  various  characteristic  plasma  spheres  of  interest.  While  much 
of  the  work  reported  here  leads  to  results  in  agreement  with  other 
authors  (Wyatt,  ^  Levine  and  Kerker^),  the  work  of  these  authors 
contains  a  considerable  number  of  errors  and  serious  ambiguities,  such 
that  it  was  felt  advisable  to  rederive  the  entire  formalism  in  detail.  An 
adequate  outline  of  the  work  and  results  contained  in  the  present  report 
is  provided  by  the  Table  of  Contents.  Part  II  will  be  concerned  with 
the  numerical  reduction  as  well  as  with  asymptotic  approximations  of 
the  analytical  results  obtained  in  Part  I. 
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MAXWELL'S  EQUATIONS  FOR  AN  INHOMOGENEOUS 


MEDIUM 


We  shall  use  MKS  units  throughout  the  present  work.  The  two  basic 
Maxwell  equations  may  then  be  written  in  the  form 


7  X  E  =  -22. 


BB 

at 


V  x  H  = 


-r 

J  +  Ht 


(1) 

(2) 


where  the  various  electromagnetic  vectors  are  related  by  the  so-called 


constitutive  relations 


B  =  |x  H 
D  =  eE 


(3) 

(4) 


and  by  Ohm's  law: 


j  =  ctE 


(5) 


? 


We  consider  the  case  where  the  magnetic  permittivity  of  the  medium  is 
that  of  the  vacuum  (p,  =  p,o),  while  the  dielectric  constant  e  and  the 
conductivity  a  may  be  functions  of  position  within  the  medium.  We 
further  assume  that  the  time  variation  of  all  electromagnetic  quantities 
has  the  form  e  1Ujt.  Substituting  this  time  variation,  as  well  as  Eq. 

(5),  into  Eqs.  (1)  and  (2),  the  two  basic  Maxwell  equations  may  be  written 
in  the  form 


7 

x  eT 

*  k2*? 

(6) 

V 

x  H 

~  -klE 

(7) 

where  k  and  k  are  given  by 
£  1 


k  =  i  p  ud 
2  o 


(8) 

(9) 


We  further  define  the  propagation  constant  k  by  means  of 


k 


2 


-klk2 


(10) 


It  is  important  to  note  that  with  the  assumptions  made  in  our  case, 
k^  is  a  constant,  while  k^  and  k  are  in  general  functions  of  position. 
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It  is  customary  to  complete  the  set  of  Maxwell  equations  with  two 
further  equations  which  give  expressions  for  the  divergences  of  B 
and  .  The  first  of  these  is 


V  •  B"  =  0 

(or  equivalently)  (11) 

?  -H  =  0 

However,  we  observe  from  Eqs.  (6)  and  (8)  that  Eq.  (11)  is  not  an  in¬ 
dependent  equation,  but  follows  automatically  from  Eq.  (6),  since  the 
divergence  of  any  curl  vanishes  identically.  The  case  of  the  divergence 
of  D  is  more  subtle  and  has  led  to  considerable  confusion  in  the  liter¬ 
ature.  In  all  cases,  we  may  write 

V  *dT  =  p  (12) 


where  p  is  the  charge  density.  The  confusion  in  much  of  the  literature 
arises  from  the  fact  that  while  it  is  true  that  we  consider  a  plasma  of 
overall  neutrality  ( "p  =  0,  i.e.  the  charge  density  of  the  free  electrons 
is  balanced  by  a  uniform  positive  background  charge,  for  example),  the 
equation  obtained  from  (12)  by  substituting  p=  0  ,  i.e. 


7  •  D  =  0 

or 

V  •  sET  =  0  ,  (13) 


which  has  been  used  by  Wyait^  and  others,  is  not  correct  for  the  case 


9 


of  time -dependent  fields.  The  correct  equation  must  be  obtained  on 
a  more  igorous  basis.  To  do  this,  we  proceed  from  the  more  funda¬ 
mental  equation  of  charge  conservation,  which  may  he  written  in  the 
form 

7  -7+  =  0  (14) 

If  we  now  substitute  Eqs.  (5)  and  (12)  and  recall  that  all  quantities  have 
the  time  variation  e  1UJt,  Eq.  (14)  becomes 

7  •  (<j  -  i<ae)E?  =0  (15) 

This  divergence  equation  is  the  correct  equation,  replacing  the 
incorrect  equation,  (13),  for  the  case  of  time -dependent  fields.  We 
note  from  Eq.  (9),  that  Eq.  (15)  can  be  rewritten  in  the  form 

V-t-k^)  =  0 

Consequently,  we  see  that  it  follows  automatically  from  Eq.  (7)  and 
thus  does  not  represent  an  independent  equation. 

Our  problem  is  thus  completely  defined  by  the  two  Maxwell  equations, 

(6)  and  (7);  the  two  divergence  equations  (11  and  15)  are  not  independent 
equations,  and  therefore  are  irrelevant. 
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III.  REDUCTION  OF  MAXWELL'S  EQUATIONS 


We  now  address  ourselves  to  the  problem  of  solving  the  Maxwell 
equations,  (6)  and  (7),  for  the  special  case  of  a  spherically  symmetric 
medium,  i.e.  one  for  which  e  and  c  (and  hence  k^  and  k  )  are 
functions  only  of  the  radial  distance  r  from  the  origin.  We  further 
assume  that  the  spherically  symmetric  medium  is  finite  in  extent, 
with  outer  radius  r  =  b. 


The  derivation  which  follows  parallels  closely  that  of  Born  and  Wolf 


(4) 


for  the  homogeneous  case;  we  also  make  use  of  the  notation  of  Wyatt. 


(1) 


Because  of  the  spherical  symmetry  of  the  problem,  it  is  most  conven¬ 
ient  to  use  spherical  coordinates  (r,  0,  <p).  Equations  (6)  and  (7)  may 
then  be  written  in  component  form  as  follows: 


k.H 
2  r 


2  . 

r  sin  0 


(rE  sin  9)  -  (rEJ 
30  cp  3<p  0 


(16) 


k2H9 


1 


ra  e 


r  sin  0 


3  «p 


£  {rVin  9) 


(17) 
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k-H 

2  cp 


dE 


£  K)  -  if 


(18) 


-k,E 
1  r 


2  . 

r  sin  9 


m 

—  (tU  sin  e\  -  ~  /rH  ) 

se  v  tp  y  atp  v  e; 


(19) 


'klE0 


1 


r  sin  0 


dH 

Up 


-  -  —  (rH  sin  0^ 
dr  \  CP  J 


(20) 


-k  E 
1  «P 


*  ("e) 


dH 

r 

"di" 


(21) 


(In  the  equations  which  follow,  we  shall  suppress  the  factor  e  1U)t, 
common  to  all  field  quantities). 


We  must  now  find  the  geneial  solution  of  the  system  of  partial  differen¬ 
tial  equations,  (16)-(21).  This  general  solution  may  be  written  as  the 

superposition  of  two  linearly  independent  field  solutions:  the  so-called 

e  -*  e  — * 

transverse  magnetic  fields  (  E  ,  H  )  for  which  =  0  ,  and  the 

transverse  electric  fields  (mE  ,  H  ),  for  which  Er  =  0  . 

Turning  our  attention  first  to  the  case  of  the  transverse  magnetic  fields 
(el5,  el3  ),  Eqs.  (16) -(21)  take  the  form 


(r  eE  sing)  -  ~  (r  ®E  ) 

96  \  <P  /  dcp  \  0/ 


=  0 


(22) 


12 


(23) 


(24) 


(25) 


I  JL 

r  dr 


(26) 


Our  problem  now  is  to  find  expressions  for  GE  ,  eH  which  satisfy  all 
of  the  Eqs .  (22)  through  (27).  This  set  of  equations  may  be  modified 
by  substituting  Eqs.  (26)  and  (27)  into  (23)  and  (24),  which  may  then  be 
written  in  the  form  (after multiplication  by  ): 
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Equation  (22)  may  be  satisfied  identically  by  choosing 


to  be  given  as  the  gradient  of  a  scalar  function 


SE  and  eE 

<p  e 


eE  =  1  ^ 

Cp  ~  r  sin  0  Sep 

e  (30) 

erp  1  a  u 

Ee  r  ae 


Moreover,  if  we  express 


in  terms  of  another  scalar  function 


as  follows: 


e 


U 


kl  Sr 


(31) 


Eq.  (3  0)  yields  the  following  expressions  for 
of  en : 


eE  and  GE 

tp  e 


in  terms 


(32) 


(33) 


Substituting  these  into  Eqs.  (26)  and  (27),  the  latter  become 


dr  59 


(34) 
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i  __a _ 

sin  0  dr  dtp 


(■  •»)  -  -  A  (■ 


6  e 

These  in  turn  may  be  satisfied  identically  if  H  and  H  are 

7  7  cp  0 

0 

expressed  in  terms  of  0  as  follows: 


i  an 

sin  0  dtp 


Of  the  original  system  of  equations,  we  have  now  satisfied  Eqs.  (22), 

(26)  and  (27),  in  the  course  of  which  we  have  obtained  expressions 

0 

(32),  (33),  (36)  and  (37)  for  the  four  angular  field  components  of  E 

0  >4  0 

and  H  in  terms  of  a  scalar  function  Q.  We  must  yet  satisfy  Eqs. 
(26)  and  (28)  through  (29)(which  are  equivalent  to  23  and  24).  To  begin, 
we  substitute  expressions  (36)  and  (37)  into  (25),  obtaining 


,  i  a  f .  „  bea\  .  i 

1  r  r  sin  0  30  \  90  J  sir 


2e 

i  a _ n 

sin  0^2 


If  we  now  substitute  Eq.  (38),  as  well  as  (36)  and  (37),  into  Eqs.  (28) 
and  (29),  these  take  the  form 
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2  2  ( 
r  sin  0  dtp  j 


(39) 


(40) 


These  equations  state  that  the  partial  derivatives  with  respect  to  tp 
and  0  of  one  and  the  same  expression  vanish  simultaneously.  Both 
equations  may  be  satisfied  simultaneously  by  assuming  that  the  bracket 
itself  vanishes,  i.e., 


+  k 


ren  + 


1 


r  sin  9  90 


1  ~2  e 

_ _ 9  0 

2  2 
r  sin  9  9? 


=  0 

(41) 


*  This  is  a  sufficient  but  not  a  necessary  condition.  The  more  general 
condition,  both  necessary  and  sufficient,  allows  the  bracket  to  be  an 
arbitrary  function  of  r.  This  in  turn  would  entail  a  modification  in 
expression  (43)  for  Er  given  below.  Here  we  shall  content  our¬ 
selves  with  the  narrower  condition  (41).  The  additional  degree  of 
freedom  provided  by  the  more  general  condition  is  somewva.t  in  the 
nature  of  a  gage  transformation. 
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_  6 
Equation  (41)  represents  an  equation  for  the  scalar  function  0.  It 

may  be  rewritten  in  the  more  concise  form 


V 


1  dk 


1  dr  v  ' 


2  e 

+  k  0  = 


(42) 


Finally,  by  substituting  Eq.  (41)  into  Eq.  (38),  the  latter  becomes  an 
6  6 

equation  for  E  in  terms  of  Q,  to  wit: 
r 


kl  5r  Sr 


r  efl 


(43) 


We  have  thus  succeeded  in  satisfying  all  of  the  original  equations 
(Eqs.  (22)  through  (27),  or  equations  derived  from  them),  in  the  course 
of  which  we  have  obtained  expressions  for  all  field  quantities  in  terms 

e 

of  a  single  scalar  function  Cl  which  must  satisfy  Eq.  (42). 


We  now  turn  our  attention  to  the  case  of  transverse  electric  fields 
(mE  ,  mH  ),  characterized  by  the  condition  mE^  s  0.  The  system 
of  equations  ( 1 6 )  through  (21)  then  takes  the  form: 


(44) 
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We  observe  that  this  system  of  equations  may  be  obtained  from  the  previous 
system  (Eqs.  (22)  through  (27))  for  the  transverse  magnetic  fields  simply 
by  the  transformations: 

E  — *  H  ,  H  — ►  -  E  ,  k1  — *  k2  ,  k2  — *  kj  (50) 

Inasmuch  as  the  transformations  in  (5  0)  convert  the  system  of  Equations 
(44)  through  (49)  identically  into  our  previous  system  of  Eqs.  (22) 
through  (27),  all  of  the  results  obtained  above  for  the  previous  system 

i 

may  be  transformed  by  means  of  Eq.  (50)  to  apply  to  the  new  system 
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of  equations  for  the  transverse  electric  fields.  Accordingly,  we  may 
define  a  scalar  function  mQ  ,  in  terms  of  which  the  transverse  electric 
field  components  are  obtained  by  transforming  Eqs.  (32),  (33),  (36), 
(37)  and  (43),  which  yields 


m 


H  =  : — 


m 


H 


m. 


H 


dr 


rn 

r  Q 


V  a 


A(^) 


k2r  sin  0  dr  dtp 


(51) 


m. 


m. 


i  d  n 


sin  0  d  tp 


dmQ 

30 


m. 


where  the  equation  satisfied  by  0  is  obtained  by  transforming  Eq. 
(42)  which  becomes 


V2  mQ  +  k2  mQ  =  0  (52) 

Note  that  in  transforming  Eqs.  (42)  and  (43),  the  term  dk^/dr  becomes 
dk^/dr  ,  which  vanishes  since  k^  is  constant. 
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0  m 

Finally,  if  we  assume  that  0  and  Q  represent  the  most  general 
solutions  of  Eqs.  (42)  and  (52),  respectively,  the  total  electromagnetic 
fields  (E  =  eE  +  mE  ,  H  =  GH  +  mH  )  in  the  medium  may  be  obtained 
by  combining  Eqs.  (32),  (33),  (36),  (37),  (43)  with  (51)  which  yields 


E  =  — 


j>L  _  _L  a  +  k* 

2  k.  dr 

dr  1  dr 


r  eQ 


(53) 


1  .2 


iTw  (r *  Sir 


(••") 


..  rm 

d  0 


(p  r  sin0  dr  dtp 


e  scp 


.  m_ 
d  Cl 

de 


(54) 


(55) 


H  ^  — 


JL_  +  k2 

^  2 
Sr 


rn 

r  Q 


H 


H 


sin 


e 

s  n 


1 _  SCQ  1  S2  /  m  \ 

n  0  dep  k  r  dr  d9  \  ) 


(56) 


(57) 


<P 


d9  k^r  sin  6  dr  dep 


(58) 


0 

The  two  scalar  functions  0  and  Q  entering  into  expressions  (53) 
through  (58)  must  be  obtained  by  solving  Eqs.  (42)  and  (52). 
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Inasmuch  as  we  have  assumed  that  depends  only  on  r  ,  both 
equations  may  be  separated  in  spherical  coordinates  according  to 
standard  methods,  and  the  general  solutions  may  be  written  in  the 
forms: 


r  en  = 


x, 

L 


[=0  m  =  0 


W(r, 


a  cosmcp  +  b  sinmcp 
m  m  T 


Pm(cos  e)  (59) 
i 


m 

r  0  = 


£  L 

m=0 


G.(r ) 


c  cosmcp  +  d  sinmcp 
m  m 


.m. 


P  (cos  e)  (60) 


where  pJn(cos  Q)  is  the  associated  Legendre  polynomial,  ^  a  ,  b  , 
i  r  mm 

c  ,  and  d  are  arbitrary  constants,  and  where  W  (r)  and  G  (r)  are 
mm  l  l 


the  general  solutions  of  the  ordinary  differential  equations 


d2W  1  dk,  dW 
_ l  _ _ l  _ l 

dr2  ^1  dr  dr 


J&U  +  1) 

2 

r 


W  =  0 

l 


(61) 


d2G 


dr 


j 1+  1 ) 

z~ 

r 


G  =  0 

l 


(62) 


Recalling  that  =  and  ^2  *S  constant»  Eq.  (62)  may  be 

rewritten  in  the  form 
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2 

K 


t.U  + 1) 


o 


(63) 


d2W 


dr 


2  _ 

k  dr 


,,  dW 
dk  l 


dr 


We  note  that  Eqs.  (42)  and  (52)  (which  lead  to  Eqs.  (61)  and  (62),  res¬ 
pectively)  were  obtained  as  a  direct  consequence  of  the  two  basic 
Maxwell  equations,  (6)  and  (7);  they  do  not  arise  from  the  require¬ 
ments  V*H  =  0  ,  V*D  =  0  ,  as  stated  by  Wyatt^);  for  example, 
the  latter  equation  is  even  incorrect,  as  we  have  seen  above. 


This  completes  the  reduction  of  Maxwell's  equations.  Given  a  specific 

dependence  k(r),  the  only  remaining  problem  is  to  find  the  general 

solutions  W  and  G  ,  and  then  to  determine  the  arbitrary  constants  in 
X  X 

expressions  (59)  and  (60)  (including  those  implicit  in  W  and  G  )  by 

Sj  l 

imposing  suitable  boundary  conditions. 
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IV.  THE  FIELD  OF  THE  INCIDENT  PLANE  WAVE 


In  the  present  report,  we  are  interested  in  the  problem  of  the  scattering 
of  a  plane  wave  by  spherically  symmetric  media  such  as  described 
above.  We  shall  assume  that  the  inhomogeneous  medium  is  finite  in 
extent,  with  outer  radius  r  =  b.  Moreover,  we  shall  assume  that  the 
incident  plane  wave  is  linearly  polarized.  The  coordinate  system  may 
then  be  chosen  such  that  the  incident  wave  propagates  in  the  positive 
z-direction  and  has  its  alectric  field  in  the  x-direction.  The  fields  of 
the  incident  wave  are  then  given  by 

i  i^z-wt) 

hi  =  e 

x 

E  1  =  E  1  =  0 

y  z 

rr  i  i(k*z- 

ri  =  -  e 

y  |A  U) 
o 

H  1  =  H  1  =  0 

x  z 
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Here  k  is  the  constant  propagation  constant  of  the  outside  medium, 
and  we  have  assumed  the  electric  field  to  be  of  unit  magnitude. 

Inasmuch  as  it  will  be  necessary  later  to  match  boundary  conditions  at 
a  spherical  surface,  we  must  reexpress  these  fields  in  spherical  coor¬ 
dinates.  For  our  purposes,  we  shall  only  require  the  r-components  of 
both  fields,  which  are  easily  shown  to  be 


ik  r  cos  ft  . 
e  sin  ft  cos  cp 


k  ik  r  cos  ft  , 

-  e  sin  ft  sin  cp 

“o® 


(65) 

(66) 


Since  it  is  much  more  convenient  to  apply  boundary  conditions  to  poten¬ 
tials  than  to  the  field  components  themselves,  it  behooves  us  to  find  the 
electric  and  magnetic  potential  functions  erj1  ,  rnQ1  from  which  the 
above  fields  may  be  derived.  To  do  this,  we  proceed  from  the  well- 
known  expansion 


I 

ik  r  cos  0  \  I  .  _  .  . 

J  =  2-  1  (2*  +  1)  J  (k  r)  P^(cos  9) 

1=1  1 


where 


Vkr)  = 


aP" 

J  2k'r 


,.i  i+7, 

2k  r  2 


J  ,1/  (kr) 


(67) 


(68) 
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Inasmuch  as 


ik^r  cos  0 


,  *  1 1 

1  ^  lk  r  cos  0 

_  be _ 

ik  r  d(cos  9) 


dP  (cos  9)  ! 

sin  9  - - -  =  P  (cos  9)  , 

d(cos  9)  * 


Eq.  (67)  may  be  rewritten  in  the  form 


ik^r  cos  0  (sin 


(2/  +  1)  j  (k*r)  P;  (cos  9) 


If  we  further  introduce  the  new  functions 


^(kIr)  =  k'rj^r)  =  4f^  JJ+l/2^r)  ,  (72) 


we  find  that  Eqs.  (65)  and  (66)  may  be  written  in  the  expanded  form 


I2  2  L _ i 

k  r  1=1 


oo 


I  l 

)  *  (k  r)  P  (cos  0) 
l  l 


sin  \ 

Hfl  =  - j-~2  )  i ^  (2j6+1)  f  (k  r)  P^(cos  9)  (74) 

V“kr 
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We  also  note  that  the  so-called  Ricatti-Bessel  functions  $  (k  r)  are 

AJ 

solutions  of  the  equation 


d  V  t2  Je(£+1) 

- £  +  k  -  _ 

,  2  2 
dr  r 


*  =  0 
’  l 


On  the  other  hand,  as  shown  in  the  preceding  section  for  a  general 

medium,  the  fields  may  also  be  derived  from  potentials  of  the  form 

(59)  and  (60)  by  means  of  the  relations  (53)  and  (56).  For  the  case  of 

constant  k  =  ,  Eqs.  (62)  and  (63)  for  G  and  W  are  both  identical 

/•  *> 

with  Eq.  (75);  consequently,  we  may  write 


°/r) 


W/r) 


r) 


Moreover,  we  may  note  from  Eqs.  (53)  and  (56)  that  only  the  m=l 
terms  in  the  general  expansions  (59)  and  (60)  are  required  to  obtain 
the  incident  fields  (73)  and  (74).  Consequently,  we  may  write  the 
electric  and  magnetic  potentials  of  the  incident  plane  wave  in  the  form 


e  i 
r  n 


=  cos  cp 


a  *  (k*r)  P  (cos  6) 


mi 

r  Q 


’  & 


b  *  (^r)  P  (cos  9) 

Xj  jo  JL 
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The  coefficients  may  be  determined  by  substituting  Eq.  (77)  into  (53) 


which  yields 


cos  cp 


a.  p,  ("os  0  ) 
l  i 


d\  + 

U  1 


(79) 


In  view  of  Eq.  (75),  this  may  be  rewritten  in  the  form 


os  cp  \ 
i  «=i 


I  i 

a  l(l+l)  4  (k  r)  P  (cos  0) 
i  it 


(80) 


By  comparing  Eq.  (80)  with  (73),  we  then  find  that  the  coefficients  a 
are  given  by 


t 


kj1  iA~1(2jR+l) 

,i2  1(1+1) 


1  ii_1(2i+l) 


k2  1(1+1) 


(81) 


such  that  the  electric  potential  (77)  becomes 


e  i 
r  n 


cos  {p  \ 

u'  L. 


V  ta1 


i£_i(2£+l )  .  „I  , 

-  4  (k  r)  P  (cos  0) 

Xj  Xj 


(82) 
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In  an  exactly  analogous  manner,  the  magnetic  potential  of  the  incident 
wave  is  found  to  be 


r 


i  sin  cp 


y-1  i 1-1  (24+1) 

4  KA+i) 


^(^r)  P,  (cos  9) 


(83) 
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V. 


SCATTERING  COEFFICIENTS  FOR  A  SPHERE 


A.  General  Case  of  an  Inhomogeneous  Sphere 

We  now  turn  to  the  general  problem  of  the  scattering  of  the  plane  wave 
described  in  Section  IV  by  a  spherically  symmetric  inhomogeneous 
medium  of  radius  i=b.  We  shall  call  the  region  outside  the  sphere 
(r  >  b)  Region  I,  that  inside  the  sphere  Region  II.  The  problem  is 
solved  by  finding  the  general  solution  for  the  fields  everywhere  and 
then  imposing  the  required  boundary  conditions. 

The  general  field  in  Region  I  is  composed  of  the  field  of  the  incident 
plane  wave  plus  the  scattered  field.  The  potentials  of  the  former  are 
given  by  Eqs.  (82)  and  (83).  The  potentials  of  the  scattered  field  are 
given  by  the  general  expression  (59)  and  (60),  where  W  (r)  and  G  (r) 
are  both  solutions  of  the  equation  (k  =  is  constant  in  Region  I) 

u,  =  0  (84) 

l 


A 


dr 


1(1  +  1) 
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Equation  (84)  has  two  linearly  independent  solutions  (one  of  these  is 

i|i  (k*r)  ,  as  we  saw  previously).  For  the  scattered  wave,  we  must 
SL 

choose  that  linear  combination  which  for  large  r  behaves  asymptoti- 
ik  r 

cally  as  e  /r  .  This  solution  is  easily  seen  to  be 


W,(r)  =  G^r)  =  HjC^Vr) 


<«> 


where  a  are  arbitr"  *y  constants. 
JL 


Moreover,  inasmuch  as  we  must  satisfy  boundary  conditions  at  r=b 
for  all  values  of  0  and  cp  ,  and  the  potentials  6Q1  ,  involve  only 

cos  cp  and  sin  cp  ,  respectively,  it  is  evident  that  the  same  must  be 
true  of  the  potentials  of  the  scattered  field  as  well  as  of  the  field  trans¬ 
mitted  into  Region  II. 


Accordingly,  the  potentials  of  the  scattered  field  in  Region  I  may  be 
conveniently  written  in  the  form 


e  s 
r  Q 


■■Zl  Y ..  ! “B1  V„  p>.  . 

k,1  A  1  ,(«+!)  4  1 

£1=1 


m  s 
r  Q 


isinCP  \  m  i ^ " 1  ( 21+ 1 )  Jl)„I  _if 
- -  )  B  -  r  (k  r)  P  (cos  0) 
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6  m 

where  the  arbitrary  constants  B  and  B.  have  been  defined  as 

9  i 

above  in  order  to  facilitate  the  later  matching  of  boundary  conditions. 
Similarly,  the  potentials  of  the  fields  in  Region  II  can  be  written  in 


m  II 
r  0 


X-l 

i  {2  HI) 
j£U+1) 


GJ.r)  P  (cos  0) 


(89) 


where  in  this  case  W  (r)  and  G  (r)  represent  particular  (rather  than 

Xj  JL 

general)  solutions  of  Eqs.  (62)  and  (63),  specifically  those  solutions 


having  no  singularities  at  the  origin. 


We  now  turn  to  the  boundary  conditions  which  must  be  satisfied  at  the 
boundary  r=b.  These  are 


(90) 


We  wish  to  express  these  boundary  conditions  in  terms  of  the  potential 

0  'JYl  _  __ 

functions  Q  and  Q.  To  do  this,  we  observe  from  Eqs.  (54),  (55), 
(57)  and  (58),  that  Eqs.  (90)  are  satisfied,  if  we  have 
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a  ,  el. 
Ta7(r  n) 


a  ,  e„II, 

“uT:  a7(r  n  1  1 

kj  (b) 


ml  m_II 

n  =  n 


,  r  =  b 


,  r  =  b 


£  <r  “nI» 


S  ,  m  II. 

a7(r  11  > 


,  r  =  b 


I  I  i  s 

where  Q  represents  the  total  potential  in  Region  I,  fl  =  fi  +  0 

Eqs.  (91)  through  (94)  represent  four  simultaneous  equations  for  the 

6  in  0  m 

four  unknown  coefficients  ,  A^  ,  B^  ,  . 


At  this  point,  it  is  convenient  to  introduce  a  new  variable  p  defined 


p  =  k  r 


Since  k  is  a  constant,  the  boundary  conditions,  Eqs.  (91)  through  (94), 
can  be  rewritten  in  the  form 


_  d  ,  el 

s7(p  n) 

ki 


_ L  ,  enn\ 

it  j.  ( p  n  )  , 

k  u(b> ap 


m  II 
Q 


p  =  x 
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p 


X 


(98) 


3  ,  mil 

—  (p  n  )  ,  P  =  x 


where  we  have  defined 


x  =  p(b)  =  kb 


(100) 


The  Ricatti-Bessel  functions  entering  into  the  potentials  (82)  and  (83) 
of  the  incident  field  and  into  the  potentials  (86)  and  (87)  of  the  scattered 
field  already  have  the  argument  p  .  We  may  also  consider  the  functions 


W  and  G  occurring  in  the  potentials  (88)  and  (89)  as  being  functions 

Xj  Xj 


of  p  ,  and  write 


e  II 

r  n 


rl 


eA  i - !2M>  W  (p)  pScos  0) 

1  Jt(jR+l)  1 


(101) 


m  II 
r  n 


l  sin 


00 

I 


■"(2—  G  (p)p1(cos  0) 
£(jt+l) 


(102) 


The  differential  equations  satisfied  by  W(p)  and  G^(p)  are  obtained 
from  Eqs.  (62)  and  (63)  by  substituting  p  =  k^r  ,  and  we  find 
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d2w 

i 


0  J  dW 
2  dn  i 


,  2  n  dp 
dp  dp 


/( i+ 1 ) 


n 


W 

i 


=  0 


(103) 


d2G 


i(X+l) 


n 


G  =  0 

i 


(104) 


where  the  relative  "refractive  index"  n  corresponding  to  a  medium 
with  propagation  constant  k  is  given  by 


n  = 


e-g.  n2  =  -j 


k 

k' 


n 


(105) 


We  now  substitute  the  potentials  (82)  through  (83),  (86)  through  (87), 

and  (101)  through  (102)  into  the  boundary  conditions  (96)  through  (99), 

0  rn 

obtaining  the  following  four  equations  for  the  coefficients  A  ,  A  , 

Xj  Xj 


m_ 

V  B 


a 


m. 


-  eA  W  (x)  =  - 

a  i 

*i(x) 

(106) 

C,(1)(x) 

Vx) 

(107) 

c,(11iX) 

eA  W>> 

*>> 

(108) 

\l 

X.kjn(b) 

34 


i 


"Bi  -  °V>>  .  -  ♦  '(,) 


where  primes  denote  differentiation  with  respect  to  the  argument. 
Multiplying  Eq.  (108)  by  k  *  and  noting  that 


(109) 


1 


k2  kl 


1 


k1II(b) 


k2  k  “(b) 


[ku(b)]i 


n02(b) 


(110) 


Eq.  (108)  may  be  rewritten  in  the  form 


*B*  ^<1)'(x> 


W  (x) 

[n2(b,]‘ 


=  *  ♦  „'(X) 


(108') 


0  xn 

Only  the  scattering  coefficients  B  ,  B  are  required  to  calculate 

x*  X 

0 

all  of  the  scattering  quantities  of  interest.  The  coefficient  B  is 
obtained  by  simultaneously  solving  Eqs.  (106)  and  (108),  which  yields 


'B 


n2  (b)  tx'(x) 


♦  A(x)  W^'(x) 


W  (x) 

l 


in  w;,x) 

<X)  W  (x) 


n22(b) 


(111) 


m 


Similarly,  B  is  obtained  by  solving  Eqs.  (107)  and  (108),  and  is 


i 


found  to  be 


35 


\ 


(112) 


Vx) 

G,(x) 


♦l(X> 


(l)i  (1)  Gt 

q  mx)  -  q }(*)  -L- 


G  (x) 

i 


The  above  expressions  may  be  written  more  compactly  by  introducing 
the  notations 


D  (x) 
l 


\  (x> 
w 


rt(x) 


0)  (x) 
l 


I 


w '(x) 
i 


W  (x) 

r 


Y/x) 


G>) 


(113) 


(114) 


(115) 


(116) 


with  which  the  scattering  coefficients  take  the  form 


e 


B 

a 


Vx) 


n22(b)  D^(x)  -  u^(x) 

(x)  -  «22(b)  r^(x) 


(117) 


36 


Yix)  -  Dj£(x) 


(118) 


tn 


tix> 


B* 


4%)  -Vx) 


This  completes  the  calculation  of  the  scattering  coefficients  for  the 


m. 


case  of  an  inhomogeneous  sphere.  (The  result  obtained  for  B  by 

X 

( 1  )  6 
Wyatt  for  the  same  case  is  incorrect.  )  The  coefficients  B  and 

JL 


m 


B  may  be  found  explicitly  for  any  given  variation  n(r)  by  solving  the 
differential  equations  (103)  and  (104).  These  were  solved  numerically 
by  Wyatt  ^  for  a  very  particular  variation  of  n(r);  they  may  be  solved 
analytically  for  n(r)  -  n^  =  constant  and  for  the  case  where  n  obeys 
a  power  law.  The  former  is  discussed  in  Section  V-B,  the  latter  in 
Section  VI-D-2  of  the  present  report. 


B.  Special  Case  of  the  Homogeneous  Sphere 

We  shall  now  specialize  the  solution  found  above  to  the  case  of  a  homo¬ 
geneous  sphere  with  a  constant  propagation  constant  ,  (k^  4  k*). 

We  thus  also  have  ^(r)  =  n^  =  constant.  This  corresponds  to  case 
(a)  in  Figure  i.  In  keeping  with  the  notation  of  Fig.  1,  the  radius  of 
the  sphere  is  now  r  =  a;  accordingly,  we  define 

y  =  kXa  (119) 
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and  the  scattering  coefficients  are  given  by  Eqs.  (Ill)  through  (112) 
with  x  replaced  by  y  . 


These  may  be  simplified,  however,  by  noting  that  W  (p)  ,  G  ( p)  for  a 

* 

homogeneous  medium  are  given  by 

W£(p)  =  G£(p)  =  yn2p)  r 


(’20) 


where  the  function  ^  i®  defined  by  Eq.  (72).  If  we  further  note  that 

Aj 


w  '(p)  =  G  '(p)  =  n2  f  (n2p) 


(121) 


where  primes  denote  differentiation  with  respect  to  the  argument,  we 
find  from  Eqs.  (113)  through  (116)  that 


<My)  =  Y/y)  =  n2  D  (n2y) 


(122) 


Accordingly,  the  scattering  coefficients  given  by  Eqs.  (117)  through 
(118)  (with  x  replaced  by  y  )  then  become 


Vy) 

n2D;.(y)  ’  Vn2y) 

D|  (n2y)  '  "2  r£(y) 

Vy) 

n2D£("2y)  -  Di(y) 

4%) 

r£(y)  '  "2D£(n2y) 

(123) 


(124) 
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These  expressions  represent  a  considerable  simplification  over  the 
corresponding  expressions  (117)  and  (118),  inasmuch  as  the  only  trans¬ 


cendental  functions  involved  in  epxressions  (123)  and  (124)  are  the  func¬ 


tions  i|r  (x)  ,  f  '(x)  ;  no  numerical  integration  of  differential  equations 

z  z 

is  required.  The  scattering  coefficients  given  by  Eqs.  (123)  and  (124) 


represent  the  complete  solution  of  case  (a)  of  Figure  (1). 
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VI  SCATTERING  COEFFICIENTS  FOR  A  SPHERE  WITH  A  CENTRAL 

CORE 


A.  General  Case 

We  now  turn  to  the  more  general  case  of  a  sphere  with  a  central  homo¬ 
geneous  core,  v.hich  includes  cases  (b),  (c)  and  (d)  of  Figure  (1).  The 
situation  to  be  considered  is  illustrated  in  Fig.  (?.)  below. 


Region  I 

FIGURE  2 


Regions  I  and  III  are  considered  to  be  homogeneous  with  constant 

,1  J  ,  III  .  ,  ...  ,11 

propagation  constants  k  and  k  ,  respectively,  while  k  may  be 
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an  arbitrary  function  of  r.  The  purpose  of  this  section  is  to  find 
expressions  for  the  scattering  coefficients  for  this  general  case;  in 
subsequent  sections,  we  shall  consider  particular  variations  of  k  * 
(or  n^  =  k^/k*  )  in  Region  II. 


As  before,  the  problem  consists  of  finding  the  potentials  everywhere 
and  then  applying  the  boundary  conditions  (96)  through  (99)  at  both 
r-a  and  r=b.  The  potentials  of  both  the  incident  and  scattered  fields 
in  Region  I  are  again  given  by  Eqs.  (82)  -  (83)  and  (86)  -  (87),  respec¬ 
tively.  Similarly,  since  k^  is  constant  in  Region  III,  the  potentials 
in  Region  III  may  be  \  ritten  in  the  form 


r 


eQm 


c  i^"1(^  +  1) 

1  £U  +  1) 


1 

Vn3  P£  ^C°S  9  ^ 


( i  r.5 ) 


m  m 
r  n 


i  sin  (p 


iH+i) 


Vn3p)  pe^cos  0) 


(126) 


Finally,  the  potentials  in  Region  II  may  be  obtained  from  the  general 
expressions  (59)  through  (60),  where  only  the  m=l  terms  can  contri¬ 
bute  because  of  the  boundary  conditions. 
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As  in  Section  V-A,  we  shall  consider  the  functions  W  and  G„  occur- 

J l  Xj 

I 

ing  in  Eqs.  (59)  and  (60)  to  be  functions  of  the  argument  p  =  k  r  ; 

W  {  p )  and  G  (  p)  then  represent  the  general  solution  of  Eqs.  (103)  and 

X'  Xj 

(104).  Inasmuch  as  Region  II  does  not  contain  the  origin,  Eqs.  (103) 
and  (104)  will  both  admit  of  two  linearly  independent  solutions,  so  that 
the  general  solutions  may  be  written  in  the  form 


W(p)  = 

YxVp)+  B;,Vp) 

(127) 

Gt(p)  - 

Y«u,<p)  +  W 

(128) 

where  •  ,  g  ,  y  .  6  are  arbitrary  constants.  The  potentials  in 

i  Xj  ij  fj 

Region  II  may  then  be  written  in  the  form 


r 


e  II 
Cl 


au+i: 


i 

P  (cos 
l 


6)[“txi(p)  + 


P.Y^P)]  (129) 


mil 
r  Q 


i  sin  <p 


I 

&=1 


l-l 


(2£+l) 


4U  +  1) 


pe(cos  e)[Y*De(p)  +  6^ve(p)]  (‘30) 


We  must  now  match  the  boundary  conditions  (96)  through  (99)  at  r-b 
(p=x)  and  r=a  (p=y).  Substituting  expressions  (82)-(83),  (86)-(b7)  and 
(129)-(130)  into  the  boundary  conditions  (96)  through  (99)  for  p  =x  , 


we  obtain  the  equations 


0 


(131) 


Vx)  +  X  ^1>(x)  -  ^Xt(x'  ■  etvx>  * 


Vx> +m^  c«1>(x) "  Yiu«(x) '  6iVx>  =  0 


(132) 


1  \  ,  e„  1 

I  ♦(I*)  +  b£  t  -  jx 

kj  kx  k  (b) 


hx>lt!(\w]=  0  <133> 


♦  >)  +mB-  C.(1)'(x)  -  Y,u'(x)  -  6.  v'(x)  =  0 


r 


'i  i'  ' 


(134) 


Similarly,  substituting  Eqs.  (129)  -  (130)  and  (125)  -  (126)  into  the 
same  boundary  conditions  for  p=v  ,  we  are  led  to  the  four  equations: 


al  X/(y)  +  "*  Yy)  =  Ct  *i(n3y) 


(135) 


Y.  U4(y)  +  V  (y)  =  d£  f^y) 


(136) 


1 


,  II. 
kj  (a) 


h  x>> +  B*  Vy)l  = 


n3  « 

c .  — rr  \lf  (n„y) 
l  ,  IL  Y  v' 


(137) 


yi'Jl{Y)  +  \  VH{V)  =  ^(n3y)  {138) 


Equations  (131)  through  (13  8)  represent  eight  simultaneous  linear 

6  m 

equations  for  the  eight  coefficients  B,  B  ,  c  ,  d  ,  <*  ,  g  , 

ij  i  tj  &  a  n 
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Y  ,  6  .  Fortunately,  these  separate  into  two  independent  sets  of 

x  x 

0 

four  equations  each.  Thus,  the  coefficient  B  may  be  found  from 

X 

Eqs.(131),  (133),  (135)  and  (137),  which  may  be  rewritten  in  the  form 


eB  CU)(x> 
i  ^X 


“i  X*W 


Yi(x) 


=  -y*> 


(139) 


x>)_  0 

\22(b)  1  n22(b) 


(140 


“1X*W  +  BiV*» 


h  x>> 1  h  Y> 


c<  W> 

n22(a)  , 

Ci  —  VV1 


(141) 


(142) 


where  Eqs.  (140)  and  (142)  were  obtained  by  multiplying  Eqs.  (133)  and 
(137)  by  and  k^(a)  ,  respectively,  and  noting  that 


h1 

k2ki! 

1 

kjk^V) 

kn  (b) 

n22(b) 

n^k^V) 

II2.  . 
n3k  (a) 

n2Z(a) 

hUI 

kb0 

2  1 

n3 

(143) 


(144) 
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Equations  (139) 
the  solution  for 


through  (142)  may  be  solved  by  standard  methods, 
0 

may  be  written  in  the  determinant  form 


and 


'B 


Vx) 


Vx) 


other  terms  as  below 


l  A, 


e.ll,w 


(i)i 

H  (x) 


-X,(x> 


l 


.Y£(x) 


—  X4(X) 


n2 


n22(b) 


YX> 


X(y) 


x/y) 


Y.(y) 


Y/y) 


-*/n3y) 


‘"z  (a>  ijfnjy) 
n„ 


The  bottom  determinant  A?  may  be  evaluated  by  first  multiplying  the 

first  row  by  -  £  ^  (x)/ £  ^  ’ (x)  and  adding  it  to  the  second  row,  which 
$/  /  & 


yields 


(145) 
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0 


(1)/ 


A2  =  C'(x) 


^U,'<x)  x’(x) 

-TT7X,(X)  -V 


c 


l 


(x)  *  n2~(b) 


C]lf,(x)  Y^(x) 

IT -  -V 


c]  '(x) 


n2  (b) 


Xi(y) 


Yy> 


W> 


x«(y) 


Yy) 


n2  (a) 


n„ 


tyn3y 


(146) 


-n  i|t  (n  y) 

If  in  the  determinant  of  Eq.  (146),  we  now  add  - —  -■  5 


times 


the  third  row  to  the  second  row,  we  obtain 


n2  (a)  tg(n3/) 


V 


•n22(a)  (1)  ■ 

~nT  h  <X,Vn3y) 


”(TT7Tx«<x) '  2 


x.(.<) 


n2  (b) 


(1)’  1 

e,  (x)  y  (x) 

—  -Y  (X)  -  4r- 


C.(1)(x) 


n2  <b> 


Xiy»- 


n3*£(n3y) 


n?  (a)t(n*y) 


■x(y) 


y  3- 


n3tin3y)  < 

Y/y)  -  T~~T - Y,(y) 


n?  <a)t(n^y) 


Y  v 


(147) 

which  may  be  easily  expanded.  Thus,  making  use  of  the  notations  (113) 
and  (114)  introduced  previously,  we  find 


46 


\ 


*2  * 


£(|1|(xl,i("jy)[r  2  I  T  2  .  1 

— t*2  (b)r*(x)X(x)  vx!]h  ,a)  wv>-  -3y^»] 

n3n2  (b)Din3y)l 

-  [n22(b)  r(x)  YJx)  -  y^(x)][n22(a)D(n3y)X^(y)  -  n3  x'(y)]  l  (148) 


The  determinant  in  the  numerator  of  Eq.  (145),  ,  is  identical  to  A^ 

(1)  (1)'  1 

except  that  £  (x)  and  £  (x)  are  replaced  by  -ijf(x)  and  -\|r  (x), 

z  z  z  z 

respectively.  Consequently,  the  value  of  A^  may  be  written  down 
immediately  from  Eq.  (148),  by  making  the  substitutions  £^(x) - ► 

Xj 


(11' 

"Vx)  ’  h  *x) 

We  then  find, 


•  i|r  (x)  ,  and  accordingly  also  T(x)  — — ♦  D  (x)  . 

z  z  z 


Vx)Vn3y)  r  2  '  ir  2  >  i 

Ai =  ~  Hn^  (b)Di(x)X;!x)  -  xHh 

n3n2  (b)IVn3y) 


•  [^2V)D(x)Yjx)  -  Y^(x)][n22(a)X  (y)D£(n3y)  -  n^y)]  I  (149) 


Thus,  the  desired  scattering  coefficient  B  becomes 
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(x)  -U(x)  -V(x) 


1 
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by  introducing  the  notations  (analogous  to  those  of  Eqs.  (J13)  through  (116): 


50 


(103)  and  (104)  may  be  solved  analytically  will  be  considered  below. 


B.  Special  Case  of  a  Spherical  Shell 


We  first  consider  the  case  of  a  spherical  shell,  for  which  the  electro¬ 
magnetic  properties  of  Regions  I  and  III  are  assumed  to  be  identical 
(k*  =  ,  and  Region  II  is  assumed  to  be  homogeneous  (n^  =  k^t*  = 


constant).  We  then  have  n^  =  k 


l/Ll .  i . 


Eqs.  (103)  and  (104)  for 


W  and  G  are  then  identical  and  have  the  previously  noted  solutions 

JC  AJ 


=  U£(  =  Vn2 p) 


(157) 


The  linearly  independent  solutions  are 


YiP)  =  V(p) 


=  Xj(n2p) 


(158) 


where  y  (z)  is  related  to  the  Neumann  function  N  L  (z)  as  follows: 
j l  1+  k 


X»  = 


V? 


w2(z> 


(159) 


The  ratios  defined  by  (154)  then  become 


^(X) 


hix) 


n2De(n2x) 


(160) 


•~f(x)  -  -  n2^  (x) 


where 


X,(x) 


(161) 
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Substituting  Eqs.  (157)  -  (158)  and  (160)  -  ( 1 6 1 )  into  expression  (155)  and  (156),  we  then  find 


I 

I 


Equations  (162)  and  (163)  then  become 


Equations  (162)  and  (163)  then  become 


(172) 


(173) 


Equations  (172)  and  (173)  -  alternatively,  Eqs.  (162)  and  (163)  -  repre¬ 
sent  the  scattering  coefficients  for  the  sphercal  shell  (case  b  of  Fig.  1). 


>ecial  Case  of  Decreasing  Refractive  Index 


We  shall  now  consider  the  case  where  the  refractive  index  in  Region  II 
decreases  with  r  ;  in  particular  we  consider  the  dependence 


(174) 


where  A  is  an  arbitrary  complex  constant.  This  corresponds  to  case 
(c)  oi  Fig.  1. 


However,  before  proceeding  to  solve  Eqs.  (103)  and  (104)  for  the  specific 
case  where  n^  is  given  by  (174),  we  shall  first  derive  some  results  of 


more  general  extent  concerning  the  analytic  nature  of  Eqs.  (103)  and 
(104).  Equation  (104)  is  of  the  general  form 


d2G 


dp 


+  f(p)G,  =  0 


(175) 


while  Eq.  (103)  has  an  additional  term  involving  the  first  derivative 


dW 


/dp  .  However,  Eq.  (103)  may  likewise  be  cast  into  the  general 


form  (175).  Thus,  if  in  (103)  we  set 


(176) 


where  p,  is  an  unknown  function  of  r  ,  we  find  with 


t  _ 


** 

II 

uw  + 

u  w 

*  4 

II 

_  II 

u  W  + 
4 

2p,'  w'  +  p,"  W 
i  it 

that  Eq.  (103)  becomes  (riter  division  by  p,  ): 


_  M 

W 

+  2 

l£  _ 

n' 

1 

-f- 

4 

n 

*  L 

nil  2n'  ,  n2 

|i  n  |i 


M  4+1) 


w  =  o 


(177) 


_  i 


The  term  proportional  to  W  may  be  eliminated  if  we  chose 

4 


Jtl 

V> 


n' 

n 


or 


p.  =  Cn 


(178) 
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I 


where  C  is  an  arbitrary  constant  which  may  be  taken  equal  to  one 
without  loss  of  generality.  Equation  (177)  then  becomes 


W 

f, 


+ 


£U+1) 

2 

P 


(179) 


which  is  indeed  of  the  form  (175),  with 


W  =  n  W 
l  t 


(180) 


Moreover,  for  the  special  case  where 


(181) 


Equation  (179)  for  W  is  identical  with  Eq.  (104)  for  G  ( p)  .  The 

Xj  Xj 

functions  n  for  which  (181)  is  satisfied  may  be  obtained  by  integrating 
Eq.  (181).  Rewriting  (181)  in  the  form 


n 

n 


ii 

r 


2 


n l 
n 


a  first  integral  is  obtained  as 

n'  =  C  n2 


C  arbitrary 


(182) 


This  in  turn  may  be  integrated  to  yield  the  general  solution 
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n 


(183) 


where  A  and  D  are  arbitrary  constants.  Thus,  for  any  n  of  the 
form  (183),  we  have 


W(p)  =  nG(p)  (184) 

We  now  return  to  the  specific  form  (174)  for  n  ,  which  is  a  particular 
case  of  (183).  Accordingly,  we  may  write 


W(p) 


h(rt 


(185) 


(since  the  equation  for  W  is  linear,  the  constant  A  in  Eq.  (184)  can 

Xj 

be  taken  as  unity  without  any  loss  of  generality),  where  G  ( p)  must 
satisfy 


d2G  [A2  -  *(£+!)] 


(186) 


The  general  solution  of  Eq.  (186)  may  be  found  by  setting 


G  =  m 
l 


(187) 


which  leads  to  the  characteristic  equation 
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m(m-l)  +  -  £U+1)J 


0 


having  the  solutions 


m 


1  ±  Vl  -  4  [A2  -  1(1+1)]* 


1,2 


Writirg 


P  =  \/  ^  +  1/2^  " 


we  may  write 


ml  =  2  +  P 


m2  =  2  ‘  P 


According  to  Eqs.  (127)  and  (128),  and  (185),  we  then  find 


u,(p> 


v,(p> 


X2(P> 


Vp) 


P+ 


-(p-^2) 


=  P 


“(P+V2) 


The  ratios  defined  by  Eq.  (154)  are  then  easily  calculated  to  be 


p  +  lk 

Vp)  =  - 

p 

p  -  lk 

Vp)  - - 

p 

p  -  V2 

M  p)  =  - 

p 

p  +  ^2 

(p)  =.  - 


> 


(196) 


Substituting  Eqs.  (192)  through  (196)  into  the  general  expressions  (155) 
through  (156),  we  find  after  some  simplifications  that  the  scattering 
coefficients  may  be  written  in  the  form 
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D.  Special  Case  of  Increasing  Refractive  Index 


1 .  Sphere  with  Central  Core 

Here  we  consider  the  case  where  the  refractive  index  n^  in  Region  II 
has  the  form 

n2  =  A  pm  (199) 


where  A  and  m  are  arbitrary  complex  constants,  with  the  exception 
that  the  case  m  =  -1  is  excluded  (this  case  was  treated  separately  in 
the  preceding  section).  We  shall  assume  that  m  is  real  and  positive 
and  that  A  has  both  positive  real  and  imaginary  parts,  such  that  (199) 
describes  an  increasing  refractive  index  (corresponding  to  curve  (d)  in 
Figure  1),  although  the  analysis  which  follows  applies  equally  well  to 
the  general  case. 


Equations  (103)  and  (104)  in  Region  II  then  take  the  form 


d2W 


,  dW 
frn  _ i 

P  dP 


a  2  2m 

A  p 


i(  i+1) 


W  =  0 


(200) 


d2G 


dp 


4* 


a  2  2  m 
A  p 


AU+1) 

2 

P 


W  =  0 


(201) 
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Both  of  these  equations  are  particular  forms  of  the  more  general 
equation 


2  r 

x  y"  +  axy'  +  (bx  +  c)y  =  0 


(202 


(51 

whose  solution  (for  r  4  0  ,  b  4  0  )  is  given  by  Kamke'  7  as 


iiiil  / 

y  =  X  2  Zv(t  VrJ) 


where 


v  ^ f  (1  -a)2  -  4c  4  0 


(203) 


(204) 


and  where  Z  is  any  one  of  the  Bessel  functions. 


v 


Thus,  Eq.  (200)  is  solved  by  making  the  identifications 


x=p,  y  =  W^,  a  =  -2m  ,  b  =  A  ,  r  =  2m  ,  c  =  -  l{fM)  (205) 


Substituting  (205)  into  (203)  and  (204)  we  accordingly  find  that  the  two 

independent  solutions  of  W  are 

Ji 
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where 


v  = 


1 

m+1 


l{JL+l)  +  (m+“2) 


(208) 


Similarly,  if  we  make  the  identifications 


x=p,  y  =  G^,  a  =  0  ,  b  =  A  ,  r  =  2m  ,  c  =  -jfcU+1)  (209) 


Eq.  (202)  becomes  Eq.  (104)  for  G  ,  and  the  two  independent  solutions 

Xj 

of  G  are  found  to  be 
l 


Ut(p) 


ip  _  /  A  m+1 

2  P 


;2io) 


where 


■  V? 


A  m+1 


(211) 


y.  = 


2  A+l 
2(m+l) 


(212) 


The  scattering  coefficients  are  obtained  by  substituting  expressions  (206) 

through  (211)  into  Eqs.  (155)  and  (156),  as  well  as  into  the  expressions 

(154)  which  enter  into  the  latter.  Inasmuch  as  the  functions  X  ,  Y  ,  U  , 

ill 

V  as  given  by  Eqs.  (206)  through  (211)  are  exceedingly  complicated, 

Xj 

no  significant  simplifications  are  possible.  We  note  also  that  for  the 


first  time  we  are  faced  with  Bessel  functions  both  whose  argument  and 
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i 


order  may  be  complex.  Moreover,  inasmuch  as  the  index  l  enters 
in  a  complicated  manner  into  the  orders  p,  and  v  of  the  Bessel 
functions,  it  is  impossible  to  write  recursion  formulas  in  t,  for  the 
Bessel  functions  involved. 

Finally,  we  note  that  the  solution  presented  in  this  section  has  already 

(2) 

been  previously  obtained  by  Levine  and  Kerker  (ci.  also  Namura  and 
Takaku^),  although  their  expressions  for  the  potentials  and  scattering 
coefficients  contain  some  errors. 


2.  Sphere  without  a  Central  Core 

Finally  we  consider  the  case  where  the  refractive  index  is  given  by 
(199)  throughout  the  entire  sphere,  i.  e.  where  the  central  core  of 
constant  refractive  index  is  absent.  This  corresponds  to  the  dotted 
curve  in  case  (d)  of  Fig.  1. 

Only  those  solutions  having  no  singularities  at  the  origin  are  then 
admissible,  and  we  have 


W(p)  =  xx(p) 

g£(p)  =  typ) 


(213) 
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where  X(p)  and  U  ( p)  are  given  by  (206)  and  (210).  The  scattering 

&  a 

coefficients  are  then  obtained  by  substituting  (213)  and  (199)  into  the 
general  expressions  (117)  and  (118),  and  we  obtain 


where 


t  ,(x) 

l 

r  a  2 ,*%  /  ,  -2m 

A  D^(x)  -  x  u^(x) 

(1), 

-2m  2 

;/  ;  x 

x  w  (x)  -  A  r  (x) 

i 

^  l  l 

v** 

"  Vx)  -  Di(x) 

•(1)(x) 

’A  '  1 

rx(x)  -  vt(x>  J 

'Vx| 

Yi(x) 


x>» 

U.'(x) 


(214) 


(215) 


(216) 


As  in  the  preceding  section,  these  expressions  are  intrinsically  quite 
complicated  and  no  further  reduction  is  possible. 
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VII.  DERIVATION  OF  THE  RADAR  CROSS-SECTION 


In  tne  preceding  sections,  we  have  obtained  exact  analytical  expressions 

6  m 

for  the  scattering  coefficient*  E  and  B  for  a  variety  of  cases. 

it  it 

Our  ultimate  interest  is  to  determine  the  radar  cross-sections  of  the 

various  objects  under  study,  and  it  is  the  aim  of  the  present  section  to 

show  that  the  radar  cross-section  can  in  general  be  expressed  entirely 

0 

in  terms  of  the  scattering  coefficients  B  and  B  . 

X  jl 

2 

The  radar  cross-section  is  defined  as  4nr  times  the  ratio  of  the 
Poynting  vector  of  the  wave  scattered  in  the  negative  z-direction  (i.  e. , 
the  radial  component  of  the  Poynting  vector  for  0=tt  )  to  the  Poynting 
vector  of  the  incident  wave  traveling  in  the  positive  z-direction.  Accord¬ 
ingly,  we  may  write  the  radar  cross-section  as 
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where  the  minus  sign  is  introduced  in  order  to  make  the  cross-section 
a  positive  quantity.  The  Poynting  vectors  entering  into  Eq.  (217)  are 


(218) 

(219) 


where  the  superscripts  i,  s,  *  stand  for  incident,  scattered  and  complex 
conjugate,  respectively,  and  where  the  subscripts  indicate  the  appropriate 
vector  component. 


We  first  turn  our  attention  to  the  Poynting  vector  of  the  scattered  wave, 
which  may  be  written  in  the  expanded  form 
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The  field  components  required  in  Eq.  (220)  are  in  general  found  from 
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expressions  (54),  (55),  (57)  and  (58)  with  Cl  and  Q  given  by  (86) 
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and  (87).  Thus,  for  example,  if  we  consider  E  ,  we  obtain 
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where  we  have  written  k  =  k  and  have  made  use  of  the  relation 
klS:  =  -k2  .  Similarly,  wa  obtain 
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We  are  generally  interested  in  the  radar  cross-section  for  the  case  where 


the  receiver  is  located  at  large  distances  from  the  scattering  object 
(far-field);  accordingly,  we  may  replace  the  Ricatti-Hankel  functions  by 
their  asymptotic  forms: 
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Furthermore,  we  are  interested  in  the  field  components  for  che  parti¬ 
cular  value  9=tt  .  Thus,  by  making  use  of  the  well-known  relations 


Substituting  these  expressions  into  Eq.  (220)  for  the  Poynting  vector, 
and  making  use  of  t.he  usual  rule  for  expressing  a  product  of  two  infinite 
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series  as  a  doubly-infinite  series,  we  obtain 


where  use  has  been  made  of  the  fact  that  =  i  |k2  |  and  where  the  bars 
denote  absolute  value. 

The  incident  Poynting  vector  II  1  is  easily  found  from  Eqs.  (219)  with 

z 

the  incident  fields  given  by  (64),  and  we  obtain 
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Substituting  expressions  (232)  and  (233)  into  (217),  we  finally  obtain  the 
desired  expression  for  the  radar  cross-section: 
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In  cor  elusion,  we  also  cite  the  corresponding  results  for  the  extinction, 
total  scattering,  and  absorption  cross-sections  (cf.  ,  for  exanple  Wyatt^): 
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VIII.  CONCLUSION 


As  described  more  fully  in  the  Introduction,  the  purpose  of  the  present 
investigation  is  to  determine  whether  the  measurement  of  radar  cross- 
section  profiles  is  a  potentially  useful  diagnostic  tool  for  ascertaining 
the  electron  density  distribution  of  inhomogeneous  plasma  spheres  of 
practical  interest.  Toward  this  end,  we  have  obtained  analytical  expres¬ 
sions  for  the  radar  cross-sections  (at  arbitrary  frequency)  of  some 
typical  examples  of  spherically  symmetric  plasma  spheres  with  increas¬ 
ing  and  decreasing  refractive  index,  as  a  function  of  radial  distance 
from  the  origin.  Specifically,  we  have  considered  the  four  different 
electron  density  distributions  illustrated  schematically  in  Figure  1. 


We  have  found  that  in  each  case,  the  radar  cross-section  is  completely 

e  m 

defined  by  means  of  two  sets  of  scattering  coefficients  B  and  B  , 

Xj  jL 

in  terms  of  which  the  radar  cross-section  can  be  calculated  by  Eq.  (234) 
of  Section  VII.  Accordingly,  we  have  obtained  analytical  expressions 


0  m 

for  the  coefficients  B  and  B  for  each  of  the  four  schematic  cases 
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illustrated  in  Fig.  1.  These  expressions  which  constitute  the  chief 
results  of  the  present  report,  are  given  by  Eqs.  (123)  and  (124)  for 
case  (a),  by  Eqs.  (172)  and  (173)  for  case  (b),  by  Eqs.  (197)  and  (198) 
for  case  (c),  and  finally  by  Eqs.  (155),  (156)  (together  with  (206)-(2 11)) 
and  (214),  (215)  for  the  discontinuous  and  continuous  distributions  of 
case  (d),  respectively. 

We  wish  to  emphasize  that  the  expressions  obtained  for  the  scattering 
coefficients  are  exact  and  are  based  on  a  full  wave  treatment  of  the 
scattering  problem,  without  recourse  to  any  mathematical  approxi¬ 
mations.  Moreover,  our  final  results  are  expressed  entirely  in  terms 
of  well-known  analytical  functions.  However,  these  expressions  -  while 
of  varying  complexity  -  are  in  general  too  cumbersome  for  hand  cal¬ 
culation.  Thus,  it  is  likely  that  a  considerable  number  of  terms  are 
required  in  the  infinite  series  (234)  defining  the  radar  cross-section 
in  order  to  achieve  the  desired  numerical  accuracy.  Moreover,  the 
problem  of  numerical  evaluation  is  complicated  by  the  fact  that  not 
all  of  the  required  Bessel  functions  are  readily  available  in  tabulated 
form.  For  these  reasons,  numerical  evaluation  of  our  results  will 
require  the  use  of  a  computer. 

This  numerical  reduction  of  our  results,  as  well  as  the  derivation  of 
asymptotic  analytical  expressions  for  the  limiting  cases  of  very  high 
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and  very  low  frequencies,  will  constitute  the  subject  matter  of  Part  II 
of  the  present  investigation. 
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